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SECTION  I 


SUMMARY  AND  CONCLUSIONS 

Worldwide  seismic  noise  levels  and  characteristics  for  1963 
are  discussed  in  this  report.  Data  for  evaluation  includes  absolute  power 
density  spectra  and  contour  maps  of  average  worldwide  microseismic  activity. 

Relative  pew  er  density  spectra  we.-e  computed  from  1963  data 
from  Worldwide  Standard  Stations.  Slopes  of  the  least-mean-square  line 
through  the  power  density  points  were  computed  and  a  pattern  of  slope  changes 
appeared  at  a  frequency  of  1.  0  cps.  A  uniform  worldwide  pattern  of  slopes 
was  observed  between  1  cps  and  2  cps.  This  suggests  two  separate  sources 
generating  microseisms  above  and  below  1  cps,  respectively,  and  that  the 
spectra  above  1  cps  are  independent  of  storms,  fronts,  etc. 

The  spectra  for  frequencies  less  than  1.  0  cps  show  greater 
seasonal  variations.  These  were  concluded  to  be  mostly  meteorological  in 
origin„ 

Monthly  contour  maps  of  average  noise  show  that  noise  is 
seasonally  variable  and  that  it  is  attenuated  at  continental  structures. 
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SECTION  II 


POWER  DENSITY  SPECTRA 


A.  PROCEDURE 

Power  density  spectra  were  computed  using  short-period 
instruments  from  selected  Worldwide  Standard  Stations.  The  data  were 
obtained  from  the  Coast  and  Geodetic  Survey  on  70  mm  film  clips.  These 
were  chosen  over  the  35  mm  film  clips  because  of  the  larger  image  size  and 
better  reproduction  quality. 

To  assure  that  the  same  method  of  determining  input  data 
was  used  throughout  the  program,  one  person  made  all  data  reductions. 

The  selected  noise  sample  was  projected  approximately  ten  times  the  original 
gram  size  on  a  large  wall-mounted  grid. 

The  availability  of  film  largely  determined  the  samples  chosen. 
Effort  was  concentrated  on  the  first  six  months  of  1963  when  it  became 
apparent  that  film  for  the  latter  months  of  the  year  would  not  be  available. 

The  noise  samples  were  objectively  chosen  to  be  representative 
of  the  recording  period.  An  attempt  was  made  to  eliminate  extremely  high 
cultural  noise  if  it  was  present  only  during  part  of  the  recording  period. 


As  far  as  was  practical,  the  horizontal  instruments  were 
sampled  at  the  same  time  that  the  corresponding  vertical  sample  was  chosen. 
Any  exceptions  were  caused  by  poor  optics  in  the  original  gram,  missing 
grams  or  reduced  quality  of  reproduction. 


Several  samples  were  taken  from  the  same  station  at  different 
times  of  the  6ame  day  to  show  the  repeatability  of  the  method.  Refer  to  the 
following  spectra  in  Appendix  A: 


Station 

Code 

Component 

Date 

Nurmijarvi,  Finland 

NUR 

N,  E 

15  January 

Quetta,  W.  Pakistan 

QUE 

Z 

24  January 

Four  samples  were  taken  on  different  days  in  the  same  month 
to  show  variations  at  a  station  during  the  month.  Refer  to  the  following: 


Station 

Valentia,  Ireland 
Kevo,  Finland 


Code 

VAL 


Component 

Z 

Z 


Date 

15  Jan.  ,  18  Jan. 

5  March,  13  March 


KEY 


Station 


Code 


Component 


Date 


Valentia,  Ireland  VAL  Z  2  April,  3  April 

Anpu,  Taiwan  ANP  Z  14  April,  28  April 

Twenty  spectra  were  computed  from  one  station  to  show  seasonal 

variations,  Valentia,  Ireland  (VAL)  was  chosen  as  the  station  because  of  film 

availability  and  the  quality  of  recordings. 

A  total  number  of  151  absolute  power  density  spectra  were 
computed  from  data  recorded  during  nine  months  in  1963. 

B.  METHOD 


The  method  used  in  the  determination  of  absolute  power  density 
spectra  was  the  technique  developed  by  Texas  Instruments  for  use  with  film 
data.  The  previously  used  polarity  technique  ^  was  improved  by  the  following: 

•  Film  traces  were  differentiated  before  infinite  clipping  to 
pre -emphasize  the  high  frequency  energy,  and 

•  Compensation  was  made  for  station  response. 

The  flow  diagram  (Figure  II- 1)  and  the  data  example  (Figure  II -2)  illustrate 
the  individual  operations  involved  in  the  method  and  their  proper  sequence. 

The  polarity  technique  applies  only  to  time  series  having  a 
Gaussian  distribution  of  amplitudes.  The  sampling  rate  on  the  I960  short- 
period  data  was  four  times  per  second  or  a  sampling  interval  of  0.  25  seconds. 
For  the  1963  study,  it  was  found  that  in  most  cases  300  seconds  of  noise  gave 
good  results  with  a  sampling  rate  of  eight  times  per  second  or  a  sampling 
interval  of  0.  125  seconds.  Figure  II  —3  compares  the  sampling  rate  and  results 
using  data  from  (a)  I960  and  (b)  1963. 

1.  Film  Measurements 

The  noise  samples  were  projected  from  film  strips  and  a 
time  sample  of  300  seconds  was  traced.  The  tracings  were  sampled  at 
intervals  of  0.  125  seconds,  yielding  2400  digital  samples  each.  A  sample 
was  recorded  as  +1  if  the  amplitude  at  that  point  was  greater  than  that  at 
the  preceding  sample  point  and  - 1  if  less.  This  corresponds  to  differentiating 
the  trace,  then  clipping  at  the  +1  and  -1  levels.  The  +l's  were  then  punched 
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1962. 


Figure  II- 1.  Flow  Diagram  for  Obtaining  Absolute  Power  Density  Spectra 
by  Polarity  Method 
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PROJECTION 


Figure  II -2.  Data  Example  for  Obtaining  Absolute  Power  Density  Spectra  by  Polarity  Method 


(a)  (b) 


Figure  XI-3.  Sampling  Rate  for  Power  Density  Spectra  Data  (a)  ,960; 

(b)  1963 

on  paper  tape  suitable  for  input  to  the  DARC  computer  for  computation  of 
relative  power  density  spectra. 

Another  set  of  measurements  was  made  on  a  selected  40- 
second  section  of  each  tracing  for  the  purpose  of  estimating  che  mean  square 
of  the  differentiated  trace.  See  paragraph  6  for  a  discussion  of  the  mean- 
square  estimate.  These  measurements  consisted  of  measuring  the  amplitude 
of  the  trace  (relative  to  any  base  line,  parallel  to  the  axis  of  the  film)  at  each 
relative  maxima  and  minima,  and  measuring  the  time  intervals  between 
successive  relative  maxima  and  minima.  For  example,  if  the  total  number 
of  maxima  and  minima  in  the  40 -second  sample  was  N,  then  N  amplitude 
and  N-l  time-interval  measurements  were  made.  These  measurements 
were  also  punched  on  paper  tape  suitable  for  computer  input. 

2.  Pre-Whitening  by  Differentiation 

Past  experience  with  the  noise  recorded  by  the  Worldwide 
Stations  has  shown  that  the  power  spectra  (before  c.llowance  for  system 
response)  generally  fall  off  rather  rapidly  with  increasing  frequency.  Since 
the  polarity  technique  gives  highest  fidelity  when  the  spectra  are  roughly 
"white,  "  the  performance  of  my  operation  tending  to  whiten  the  noise  before 
infinite  clipping  would  improve  the  spectral  estimates.  One  such  operation 
is  differentiation,  which  corresponds  to  a  6  db/ octave  multiplication  of  the 
power  spectrum.  With  digitiz'd  data,  differentiation  maybe  approximated 
by  subtracting  the  previous  sample  from  a  given  sample.  Symbolically, 
the  operation  is 

At  g'(t)  =  g(t)  -  g(t  -  At) 

This  operation  is  equivalent  to  convolution  with  a  two-point  (+1,  -1)  operator 
whose  power  response  is  2{i  -  cos  2nf At) . 


n- 


3.  Computation  and  Correction  of  Polarity  Autocorrelations 


The  polarity  autocorrelations  were  computed  as 
2400-k 

$  {kAt)  =  E  h(nAt)  h(nAt  +  kAt)  for  k  =  0,  1,2, . .  250 

n=l 

where  h(n/,t)  is  either  +1  or  -1,  and  At  =  0.  125  seco.id. 


The  autocorrelations  were  then  "arc-sine  corrected"  to 
obtain  estimates  of  the  autocorrelations  that  would  have  been  obtained  if  the 
data  had  been  fully  quantized.  This  correction  is 


$'(kAt)  =  sin  f  for  k  =  0,  1,  2, . .  250. 

f 

The  theory  of  polarity  autocorrelations  and  the  arc-sine 
correction  has  been  fully  presented  in  a  previous  reporc2  and  thus  will  not 
be  repeated  here. 

4.  Computation  of  Relative  Power  Density  Spectra 


The  corrected  polarity  autocorrelations  were  Fourier  trans¬ 
formed  as  follows: 

250 


$(jAf)  =  2  $'{kAt)  cos  2nvjAf)  (kAt) 

k  -  250  1 


At 


j  =  0,  1,  2, . .  125,  where  Af  = 


1 


500At  ’ 


These  Fourier  transforms  were  then  Harmed  by  convolving  with  a  three-point 
(1/4,  l/2,  1/4)  operator  to  obtain  the  smoothed  power  density  spectra.  These 
spectra  are  relative  spectra  only,  since  the  §'(kAt)  are  normalized  autocorrela¬ 
tion  functions. 


5.  Post- Whitening  and  Compensation  fo”  System  Response 

These  two  operations  were  performed  simultaneously  for 
reasons  of  computer  efficiency  and  precision.  Since  the  power  response  of 
the  whitening  operation  is  2(1  -  cos  2nfAt),  post -whitening  to  remote  this 
effect  consists  of  dividing  the  spectrum  by  2(1  -  cos  2rrfAt).  Compensation 
for  system  response  consists  of  dividing  the  spectra  by  H2(f)  where  H(f)  is 
the  amplitude  response  of  the  system.  Amplitude  response  was  obtained 
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from  U.  S.  Department  of  Commerce  bulletin  "Instrumentation  of  the  World- 
Wide  Seismograph  System,  "  Figure  II -4.  This  system  response  includes 
all  effects  of  the  recoiding  apparatus  (seismometers^  galvanometers,  etc.) 
and  relates  absolute  ground  motion  to  film  trace  deflection. 

6.  Mean  Square  Estimation 

To  obtain  absolute  power  density  spectra  from  the  relative 
power  density  spectra  obtained  with  normalized  autocorrelations,  it  is  neces¬ 
sary  to  multiply  each  by  the  mean  square  of  the  corresponding  time  function. 
Using  polarity  data,  we  have  no  way  of  computing  the  mean  square;  however, 
we  may  obtain  a  good  estimate  of  the  mean  square  of  each  whitened  trace.  We 
approximate  the  unwhitened  film  trace  between  successive  relative  maxima 
and  minima  by  half  cycles  of  cosine  functions  of  appropriate  amplitudes  and 
frequencies.  This  approximation  is  illustrated  by  the  g.(t)  in  the  diagram 
below:  ,  1 


A.  +  A  A.  -  A 

Then  g.(t)  =  — - - -  +  — - - -  cos  (n  )  for  the  interval  6t. . 

x  c  c  6t.  l 

l 


Differentiating  g.(t)  we  obtain 


i 

Ai  +  1  '  Ai  TT  t 

g'.(t)  = - - -  sin  (tt  — )  for  the  interval  6t.. 

i  c.  ot.  ot.  1 

l  l 

Then  the  mean  square  of  g'(t)  may  be  estimated  as 

fit. 


E-M-s-=rWf  / 

i  i  1  t  =  o 

A  A  X2  2 

A,  -  A.)  TT 

r—  \ ,  c*  1  r*  (  1  +  1  1 

E.  M.  S.  =  —  — -  E - 


dt 


E  fit.  i 
i  1 


46t. 


6ti  2  t 

/  ‘8in  <n6T)dt 

t  =  o  i 


.  n2(A.  -  A.)‘ 

E.M.S.  =  =rrT“  E  —  1 


E  fit.  i 

i  * 


86t. 

i 
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MAGNIFICATION 


Figure  II-4.  Frequency  Response  of  the  USC&GS  World-Wide  Standard 
Short  Period  Seismograph  System 


The  and  6t^  measurements  were  made  directly  from  the  film  as  described 

in  paragraph  1.  These  measurements  were  input  to  DARC,  and  estimated 
mean  squares  for  each  whitened  noise  sample  were  computed  using  the 
formula  just  derived. 

7.  Conversion  of  Relative  Spectra  to  Absolute  Spectra 

Absolute  power  density  spectra  were  obtained  by  multiplication 
of  the  relative  power  density  spectra,  previously  corrected  for  whitening 
and  system  response,  by  their  respective  mean-square  estimates,  and  by  a 
constant.  This  constant  contains  such  scaling  factors  as: 

•  Magnification  factor  involved  in  projecting  film  and 
measuring  A.'s 

•  Variations  in  gain  from  that  gain  at  which  the  system 
response  was  calibrated,  and 

•  Computer  scale  factors  introduced  in  data  processing. 

8.  Presentation  of  Spectra 

The  logarithm  of  each  point  of  each  absolute  power  density 
spectrum  w:s  written  on  magnetic  tape  by  the  DARC  computer  in  a  special 
configuration.  This  configuration  was  such  as  to  provide  a  plot  of  the 
logarithm  of  the  power  density  function  versus  linear  frequency  by  playback 
of  the  magnetic  tape  through  a  digital-to -analog  converter  and  oscillograph. 


SECTION  III 


EVALUATION  OF  POWER  DENSITY  SPECTRA 

The  power  density  spectra,  computed  on  a  worldwide  basis 
and  converted  to  absolute  values,  lends  itself  to  a  study  of  the  noise  of  the 
world. 


The  spectra  plots  are  located  in  Appendix  A.  In  each  case, 
the  slope  of  the  plots  were  determined  by  the  least  mean  square  method. 

It  is  interesting  to  note  the  change  in  siope  on  the  plots  in  the  range  of  1.  0 
cps.  This  slope  change  is  interpreted  to  suggest  that  two  source  mechanisms 
are  causing  the  noise  above  and  below  1.  0  cps.  The  same  change  in  slope  was 
observed  by  Vinnik  and  Pruchkina  (1964). 

Several  comparisons  were  made  between  the  slopes  obtained 
from  the  spectra  plots.  The  first  is  the  number  of  stations  whose  spectral 
slope  (frequency  1.  0  cps  or  greater)  falls  within  each  slope  increment  of 
10  db/octave  (Figure  III  - 1 ) .  It  was  noted  that  84.  1  percent  of  the  slopes' 
values  fall  within  the  range  of  10-40  db/octave  and  that  of  these  48.  8  percent 
fall  within  the  range  of  20-30  db/octave.  This  worldwide  uniformity  in  noise 
betweer  the  periods  of  0.  5  sec  to  1.  0  sec  suggests  a  constant  universal  noise. 

Several  theories  have  been  advanced  thus  far  as  to  the  reason 
for  the  presence  of  microseisms  throughout  the  world.  Iyer  (1962)  suggests 
that  the  earth  itself  is  filled  with  noise.  Other  causes  from  cultural  noise 
to  sea  storms  have  been  investigated.  The  problem  lends  itself  for  f\  rther 
study. 


The  power  density  spectra,  plots  yield  an  increase  in  slope 
between  0.  25  c^s  to  1.  0  cps.  This  spectral  range  corresponds  with  the  theory 
that  microseisms  in  the  2-6  sec  period  range  receive  their  energy  Trom 
pressure  fronts. 

A  comparison  was  made  between  the  slopes  of  the  spectra 
plots  and  the  station  distance  from  large  bodies  of  water.  Station  geologic 
structures  were  also  taken  into  consideration  (refer  to  Figures  III -2  and  III-3). 
As  can  be  seen  from  Figure  III -2 ,  the  noise  present  in  the  3  sec  period 
range  levels  off  at  a  distance  of  around  300  miles  from  the  shore.  The  high 
frequency  spectra  (Figure  III -3)  drops  off  much  more  rapidly  although  it 
apparently  levels  off  at  about  300  mi’'s  also.  There  is  no  apparent  relation¬ 
ship  cf  these  noise  spectra  with  the  geologic  environment  near  the  recording 
stations. 


The  results  from  the  seasonal  study  at  Valentda,  Ireland  show 
that  the  seasonal  change  in  the  spectra  around  1.  0  cps  and  above  is  less  than 
that  below  1.  0  cps. 
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SECTION  IV 


VISUAL  NOISE  MEASUREMENTS 

In  addition  to  the  noise  power  density  spectra  study,  measure¬ 
ments  and  companions  of  average  microseismic  levels  on  a  worldwide  basis 
were  made.  Visual  noise  measurements  were  made  at  selected  stations 
in  order  that  noise  levels  could  be  compared.  To  improve  station  distribution, 
35  mm  records  from  the  Canadian  Network  were  obtained  and  used  along  with 
Worldwide  Standard  Station  records. 

Measurements  of  period  and  amplitude  of  the  average  maximum 
short-period  and  long-period  noise  on  the  respective  instruments  were  made 
at  each  station.  A  comparison  was  made  between  one  measurement  each  day, 
one  measurement  every  other  day,  one  measurement  every  third  day  and 
one  measurement  every  fifth  day  at  several  stations.  It  was  decided  that 
ten  measurements  each  month  could  be  made  and  not  significantly  affect  the 
average  maximum  measurements. 

The  measurements  were  taken  at  the  same  time  at  each  station, 
within  limits.  Anomalous  data,  e.  g.  ,  long-period  noise  on  short-period 
instruments,  were  eliminated  by  limiting  the  measurements  to  noise  in  the 
0.  5  to  2.  0  sec  period  range  as  the  short-period  instruments  and  the  3.  0 
to  8.  0  sec  period  range  on  the  long-period  instruments. 

The  measurements  were  converted  to  ground  motion  using 
the  approximate  response  curves  (Figure  IV- 1  -  USC&GS  Standard  Stations; 
Figure  IV-2  -  Canadian  Network  Stations)  and  averaged  for  each  month. 

Contour  map?  of  the  average  maximum  noise  for  each  month 
were  constructed.  These  maps  are  contained  in  Appendix  B.  Only  nine 
months  of  the  year  were  included  as  data  for  July,  August  and  September 
were  not  available  when  data  reduction  was  terminated. 

To  assure  consisted  i  measuring  techniques,  the  measure¬ 
ments  were  made  by  one  person,  and  one  analyst  was  responsible  for  data 
reduction  and  the  preparing  of  the  maps. 

Several  large  areas  of  the  world  have  sparse  data  available 
for  this  study.  An  attempt  has  been  made  to  subjectively  qualify  the  maps. 
Solid  lines  indicate  adequate  control  points  and  dashed  lines  indicate  a  lack 
of  control  points. 

As  can  be  seen  from  the  maps,  the  noise  level  increases  during 
the  winter  months  and,  as  was  expected,  the  noise  level  was  attenuated  at 
continental  structures. 
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Figure  IV- 1.  Frequency  Response  of  the  USC&GS  World-Wide  Standard 
Seismograph  Systems 
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Figure  IV-2,  Frequency  Response  of  the  Canadian  Network  Seismograph 
Systems  {Vertical  Components) 


Microseisms  studies  are  needed  in  conjunction  with  seismicity- 
studies  to  obtain  a  better  understanding  of  the  noise  that  can  be  expected  at 
various  stations  around  the  world.  This  noise  can  be  used  to  calculate  the 
theoretical  limits  of  perceptibility  for  e?ch  station. 

Microseisms  studies  should  also  add  to  the  knowledge  of 
seasonal  variations  at  the  stations  studied.  This  information  could  aid  in 
determining  what  stations,  if  any,  should  change  their  gains  to  cope  with 
these  variations. 


The  1963  microseisms  study  yielded  the  following  results 
as  to  the  gain  and  noise  at  the  stations  studied. 

List  of  Stations  with  Cultural  Noise  Limiting  the  Gain 


Addis  Ababa,  Ethiopia 
Ann  Arbor,  Michigan 
Athens,  Greece 
Baguio,  Philippines 
Bogota,  Colombia 
Helwan,  Egypt 
Hong  Kong 


La  Paz,  Bolivia 
Lubbock,  Texas 
Malaga,  Spain 
Minneapolis,  Minnesota 
Quito,  Ecuador 
Rabaul,  New  Bri„  tin 
Stuttgart,  West  Germany 


List  of  Stations  with  Noise  in  the  4-8 

Albuquerque,  New  Mexico 
Blacksburg,  Virginia 
Bulawayo,  Southern  Rhodesia 
Chiengmai,  Thailand 
Golden,  Colorado 
Kipapa,  Hawaii 

List  of  Stations  Where  C 


Second  Range  on  SP  Instruments 

Kongsberg,  Norway 
State  College,  Pennsylvania 
South  Pole,  Antarctica 
Tasmanian  University,  Tasmania 
Tucson,  Arizona 

ains  Could  Be  F.aised 


Chiengmai,  Thailand  (LP  only) 

Hong  Kong 

Nairobi,  Kenya  (only  if  better  photographic  technique  employed) 

New  Delhi,  India 

Any  further  work  in  this  area  should  include  more  stations  as 
they  become  available  and,  if  no  World  Standard  Stations  are  a  ailable  for 
areas  such  as  Eastern  South  America  and  Western  Africa,  reliably  calibrated 
non-Standard  station  data  should  be  used. 

The  overall  quality  of  the  film  reproductions  used  in  this  study 
varied  widely  from  excellent  to,  in  a  few  cases,  unuseable.  The  reasons 
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f  it  some  records  were  unuseable  could,  in  some  cases,  be  attributed  to 
j.  -  reproductions  while  others  were  due  to  poor  original  grams.  In  addition, 
n.  v  stations  had  excessive  trace  width  which  made  accurate  measurements 
of  low  noise  levels  difficult. 

The  most  common  problem  of  the  Wirld  Standa:  J  Network 
records,  as  a  whole,  is  excessive  drift  on  the  long-period  horizontal  instru¬ 
ments,  possibly  caused  by  inadequately  insulated  vaults. 
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SECTION  V 


RECOMMENDATIONS 

1.  Noise  background  studies  are  greatly  improved  if  a  technique  of 
obtaining  power  density  spectra  is  used. 

2.  Continued  '"orldwide  effort  should  include  gathering  of  meteorologies 
information  to  aid  in  the  interpretation  of  spectra. 

3.  Any  additional  effort  to  investigate  worldwide  microseisms  should 
have  as  a  prime  objective  the  association  of  spectra  and  observations  with 
existing  theories  of  microseisms  origin. 

4.  Additional  theoretical  and  practical  work  should  be  associated  to 
explain  constant  microseismic  background  apparently  not  associated  with 
storms,  fronts,  etc. 
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ABSOLUTE  POWER  DENSITY  SPECTRA 


APPENDIX  A 


ABSOLUTE  POWER  DENSITY  SPECTRA 

This  appendix  contains  absolute  power  density  spectra  obtained 
from  the  1963  data.  The  station  abbreviations  are  listed  with  the  location  of 
each  station  in  Table  A-l.  The  spectra  are  presented  in  Figures  A- la  through 
A- If.  The  Hat  of  at&tiona  used  and  the  dates  of  the  records  are  presented  in 
Tables  A-2a  through  A-2f  for  the  associated  spectr?.!  illustrations. 


TABLE  A-l 


STATION  ABBREVIATIONS 


STATION 

LOCAUON  STATION 

LOCATION 

AAE 

Addis  Ababa,  Ethiopia 

KIP 

Kipapa,  Hawaii 

ADE 

Adelaide,  S.  Australia 

KON. 

Kongsberg,  Norway 

AFI 

Afiamalu,  W.  Samoa 

LON 

Longmire,  Washington 

ALQ 

Albuquerque,  New  Mexico 

MAL 

Malaga,  Spain 

ANP 

Anpu,  Taiwan 

MAN 

Manila,  Philippines 

ATU 

Athens,  Greece 

MDS 

Madison,  Wisconsin 

BAG 

Baguio,  Philippines 

MUN 

Mundaring,  W.  Australia 

BKS 

Berkeley,  California 

NAI 

Nairobi,  Kenya 

BLA 

Blacksburg,  Virginia 

NUR 

Nurmi jarvi,  Finland 

BUL 

Bulawayo,  5.  Rhodesia 

PLM 

Palomar,  California 

CCG 

Camp  Century,  Greenland 

PMG 

Port  Mi  resby,  New  Guinea 

CHG 

Chiengmai,  Thailand 

PRE 

Pretoria,  S.  Africa 

CMC 

Copper  Mine,  Canada 

PTO 

Porto,  Portugal 

COP 

Copenhagen,  Denmark 

QUE 

Quetta,  W.  Pakistan 

COR 

Corvallis,  Oregon 

SCP 

State  College,  Pennsylvania 

GDH 

Godhavn,  Greenland 

SEO 

Seoul,  S.  Korea 

GOL 

Golden,  Colorado 

SHL 

Shillong,  India 

GSC 

Goldstone,  California 

SPA 

South  Pole,  Antarctica 

GUA 

Guam,  Mariana  Islands 

TOL 

Toledo,  Spain 

HUR 

Honiara,  Solomon  Islands 

VAL 

Valentia,  Ireland 

1ST 

Istanbul,  Turkey 

WES 

Weston,  Massachusetts 

KEV 

Kevo,  Finland 

WIN 

Windhoek,  S.  Africa 

TABLE  A-2a 


ABSOLUTE  POWER  DENSITY  SPECTRA  LOCATED  IN  FIGURE  A-la 


STATION 

DATE 

COMPONENT 

GAIN(K) 

AAE 

16  January 

SPZ 

25 

ADE 

15  January 

SPZ 

25 

ADE 

15  January 

SPN 

25 

ADE 

15  January 

SPE 

2_> 

ADE 

16  April 

SPZ 

25 

AFI 

16  January 

SPZ 

12.  5 

ALQ 

7  January 

SPZ 

200 

ALQ 

19  April 

SPZ 

400 

ALQ 

19  April 

SPN 

400 

ALQ 

19  April 

SPE 

400 

ANP 

14  April 

SPZ 

6.  25 

ANP 

28  April 

SPZ 

6.  25 

ANP 

28  April 

SPN 

6.  25 

ANP 

28  April 

SPE 

6.  25 

ATU 

16  January 

SPZ 

12.  5 

A.U 

19  April 

SPZ 

12.  5 

ATU 

19  April 

SPN 

12.  5 

ATU 

19  April 

SPE 

12.  5 

BAG 

21  January 

SPZ 

25 

BAG 

21  January 

SPN 

25 

BAG 

21  January 

SPE 

25 

BAG 

19  April 

SPZ 

25 

BAG 

19  April 

SPN 

25 

BAG 

19  April 

SPE 

25 

BKS 

10  January 

SPZ 

25 
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Figure  A- la.  Absolute  Power  Density  Spectra  Ob 
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TABLE  A-2b 


ABSOLUTE  POWER  DENSITY  SPECTRA  LOCATED  IN  FIGURE  A- lb 


STATION 

DATE 

COMPONENT 

GAIN(K) 

BKS 

19  April 

SPZ 

25 

BKS 

19  April 

SPN 

25 

BKS 

19  April 

SPE 

25 

BLA 

14  January 

SPZ 

50 

BLA 

10  April 

SPZ 

25 

BUL 

.  15  April 

SPZ 

100 

CCG 

7  April 

SPZ 

100 

CHG 

18  March 

SPZ 

400 

CHG 

7  June 

SPZ 

400 

CHG 

7  June 

SPN 

400 

CHG 

7  June 

SPE 

400 

CMC 

1  May 

SPZ 

ZOO 

CMC 

1  May 

SPN 

ZOO 

CMC 

1  May 

SPE 

ZOO 

COP 

15  January 

SPZ 

12.  5 

COP 

19  April 

SPZ 

12.  5 

COR 

5  January 

SPZ 

25 

COR 

1  April 

SPZ 

25 

GDH 

16  January 

SPZ 

12.  5 

GOL 

16  January 

SPZ 

400 

GOL 

5  April 

SPZ 

400 

GSC 

3  April 

SPZ 

200 

GUA 

26  May 
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6.  25 

GUA 

26  May 

SPN 

6.25 

GUA 

26  May 

*  SPE 

6.25 

HNR 

16  January 

SPZ 

12.  5 

HNR 

16  January 

SPN 

12.  5 

r 

f 


j i 
s 


1 


I 

I 

I 

f 


r 


M.  %. 


I 

i 

I 

! 

! 

I 

1 

I 

i 

! 


a 

u 


1C 


10' 


Frequency  in  cpa 

BKS  19  APR, 

SPZ  25  K 


c 

1C  - 


2 

1 


Frequency  in  cpa 


19  APR. 

SPN  25  K 


BXS  19  APR. 

SPE  25  K 


V 


spz  50  K 


spz  25  K 


m 

a 

v 

n 

u 

V  m 


1 0 


10 


10 


A 


-47.  41 


V 


\ 


o  r  ^ 

a,  £  10 


a 

o  10 

4) 

.O 


\\  \ 


'  ^ 


\, 


-25.  87 


U, 

V 


V’v‘ 


.25  .50 

»L 

SP/, 


1.0 


— i _ 

1.  5 


1  requency  in  cpe 


~z!o 


1  5  APR. 
100  K 


IS  APR.  CHG  7  JUNE  CMC 

100  K  SPZ  400  K  SPZ 


b 


Absolute  Power  Density  Absolute  Power  Density  Absolute  Power  Density 


-68. 70 


X  0 

[ho  - 


-26.  98 


.  25  .  50  1.0  1.5 

Frequency  in  cpa 

CHG  7  JUNE 

SPN  400  K 


10 


c  1 

X  10  - 

Q 

V  0 

ry°  - 
»■ E 

«  10  - 
3 

0  -> 
.2  io- 


■54. 78 


.  25  .  50  1.  0  1.5 

Frequency  in  cpa 

CMC  1  MAY 

SPN  200  K 


-39.  41 


•57. 29 


.25  .50  1.0  1.5 

Frequency  in  cpa 


7  JUNE 
400  K 


.  25  .  50  1.  0  1.5 

Frequency  m  cpa 


I  MAY 
200  K 


52.95 


.  25  .50  1.0  1.5 

Frequency  in  cpa 


1  MAY 
200  K 


v  ■> 

Q  io2 

U 

*rga.  1 

&  E10 


3  10  - 

0 

.2  -1 

<  10  - 


•h 


■31.  07 


1  \\,  *19-J6 


* 


25  .  50  1.0  1.5 

Frequency  in  cpa 


15  JAN. 
12.  5  K 


1C  - 

3 

« 

10  - 

c 

0) 

Q 

2 

Ih 

10  - 

K 

O  c 

a,  E 

io1  - 

4‘ 

3 

0 

rH 

0 

10  - 

CO 

< 

1  o" 1  - 

\ 


-61.  32 


-18.  02 


\ 


.  25 

COP 

SPZ 


-J _ 1 - 1 - I - -I 

50  1.0  1.  5  2.0 

Frequency  in  cps 

19  APR. 

12.  5  K 


10 


6 


I  !05  - 


\ 

-4  •*  11 

•  \ 

m'-t\ 

\\  ■25•8, 
'Av 

'  % 

1'  ■ 

1 

.25  .50 

1.0  1.5 

Frequency  i_  cps 

GDH 

16  JAN, 

SPZ 

12.5  K 

2.0 


c 

V 

a 


V 

3 

O 

» 

< 


,4  -  ^ 

v\ 


SPZ  15  K 


10 


10 


H  10 

c 

v 
Q 

u 

0 

£  E  10 


3 

r-4 

O 

n 

.O 


\  w 


-77. 45 


.2  10 


-1 


io'2- 


Frequency  in  cps 

GOL  16  JAN. 

SPZ  400  K 


m 

a 

V 

Q 


4 

10 


3 

10 


j!  2 

•tV0 
0  ? 
a,  c 

2  10 


o 

« 

-O 

« 


10 


COR  1  APR. 

SPZ  25  K 


a 

C 

v 

Q 


v 

«-» 

3 

O 

u 

.C 

< 


GOL  5  APR. 

SPZ  400  K 


4 


SPZ 


10 


a  .  .4 

D  10  - 

Q 

u  3 

*v°  • 

£  £  , 

2 

10  - 


V 


V 

4-» 

3 
. — . 
o 

c 


10 


.25  .  50  1.0  1.5 

Frequency  in  cpa 


2.0 


200  K 


GUA 

SPE 


26  MAY 
6.25  K 


a 

a 

V 

Q 


10 


10 


10‘ 


u 
m 
>N 

S,  lio* 


i 


~  10  - 


10“ 


A 


-160.  40 


.25  .50  1.0  1.5 

Frequency  in  cpa 

GUA 
SPZ 


2.0 


26  MAY 
6.25  K 


>»  10 


c 

v 

Q 

u 


10 


£  alO 
0  c 
o.  c 


2  10 


V 

0 

■ 

45 

< 


10 


-28. 87 


F r.-quency  In  cpa 


HNR 

SPZ 


16  JAN. 
12.  5  K 


*  10  - 


C  4 
0  10 

>V°3 

v  .  -2 

M 

3 
0 


4) 

< 


10 

101 


Frequency  in  cpa 


GUA 

SPN 


10 


a 

V 

a 

k 


0  F' 

0,  C 

V 

4-* 

P 

0 

a 


X 


26  MAY 

6.25K  SPN  12. 5  K 

Figure  A- lb.  Absolute  Power  Density  Spectra  Obtained  From  1963  Data 
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TABLE  A-2c 


ABSOLUTE  POWER  DENSITY  SPECTRA  LOCATED  IN  FIGURE  A-lc 


STATION  DATE  COMPONENT  GAIN(K) 


HNR 

4  March 

SPZ 

12.  5 

1ST 

2  January 

SPZ 

25 

1ST 

19  January 

SPZ 

25 

1ST 

20  April 

SPZ 

25 

1ST 

20  April 

SPN 

25 

1ST 

20  April 

SPE 

25 

KEV 

12  January 

SPZ 

25 

KEV 

5  March 

SPZ 

25 

KEV 

13  March 

SPZ 

25 

KIP 
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SPZ 
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KIP 
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50 
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Figure  A-lc.  Absolute  Power  Density  Spectra  Obtained  From  1963  Data 


1  ABLE  A-2d 


ABSOLUTE  POWER  DENSITY  SPECTRA  LOCATED  IN  FIGURE  A- Id 


STATION  DATE  COMPONENT  GAIN(K) 


MAN 

3  June 

SPZ 

12.  5 

MAN 

3  June 

SPN 

12.  5 

MAN 

3  June 

SPZ 

12.  5 

MDS 

18  January 

SPZ 

100 

MUN 

18  January 

SPZ 

25 

MUN 

18  January 

SPN 
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MUN 

18  January 

SPE 

25 

MUN 

5  April 
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25 

NAI 
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APPENDIX  B 


NOISE  MAPS 

This  appendix  contains  the  noise  maps  for  nine  months  in  1963. 
The  short  period  maps  (0.  5  -  2.  0  sec)  precede  the  long  period  maps  (3. 

8.  0  sec)  for  each  month  so  that  comparison  may  be  more  easily  made, 
maps  comprise  Figures  B-l  through  B-18. 
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Figure  B-2.  World  Map  of  3.  0-8.  0  Second  Microseismic  Activity,  January,  1963 
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Figure  B-3.  World  Map 


j  Figure  B-4.  World  Map  of  3.  0-8.  0  Second  Microseismic  Activity,  February,  1963 


Figure  B-5.  World  Maj 
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Figure  B-6.  World  Map  of  3,  0-8.  0  Second  Microseismic  Activity,  March,  1963 
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Figure  B-7.  World  Map  c 


B-8 


Figure  B-9.  World  Ma 
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Figure  B-9.  World  Map  of  0.  5-2.  0  Second  Microseismic  Activity,  May,  1963 
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Figure  B-10.  World  Map  of  3.  0-8.  0  Second  Microseismic  Activity,  May,  1963 
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Figure  B-ll.  World  Map  of  ( 


Figure  B-ll.  World  Map  of  0.  5-2.  0  Second  Microseismic  Activity,  June,  1963 
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Figure  B-12.  World  Map  of  3.  0~8.  0  Second  Microseismic  Activity,  June,  1963 
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B-13 


Figure  B-13 


^  Figure  B- 14.  World  Map  of  3.  0-8.  0  Second  Microseismic  Activity,  October,  1963  | 
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Figure  B-  16.  V 


Figure  B-17.  World  Map  of  0.  5] 
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